1 



(N 
O 
O 
(N 



(N 



> 

in 
m 

o 

(N 
O 

o 



The Nature of Dark Matter in Elliptical Galaxies: Chandra Observations 
of NGC 4636 
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Code 662, Greenbelt, MD 20771, USA 

We determine the total enclosed mass profile from 0.7 to 35 kpc in the elliptical galaxy NGC 4636 based on 
the hot interstellar medium temperature profile measured using the Chandra X-ray Observatory, and other X-ray 
and optical data. The total mass increases as r^'^ to a good approximation over this range in radii, attaining a 
total of ~ 1.5 X 10^^ Mq (corresponding to Mtot/Lv ~ 40) at 35 kpc. We find that at least half, and as much 
as 80%, of the mass within the optical half-light radius is non-luminous, implying an exceptionally low baryon 
fraction in NGC 4636. The large inferred dark matter concentration and central dark matter density, consistent 
with the upper end of the range expected for standard cold dark matter halos, imply that mechanisms proposed 
to explain low dark matter densities in less massive galaxies are not effective in elliptical galaxies. 



1. Context 

According to recent estimates, 80-90% of the 
matter in the universe is non-baryonic. With 
the presence of extended dark matter halos in 
galaxies of all morphological types now well- 
established, attention is focusing on comparing 
the detailed mass distribution with theoretical 
predictions of galactic dark halo structure. The 
shape of the dark matter distribution is deter- 
mined by the initial density perturbation spec- 
trum, the coupled dynamical evolution of bary- 
onic and non-baryonic constituents, and the na- 
ture of the dark matter itself. Therefore, its mea- 
surement represents a powerful diagnostic of fun- 
damental astrophysical processes and parameters. 

The standard cold dark matter (CDM) model 
is highly successful in explaining the distribution 
of mass in the universe on scales ranging from 
galaxies on up, but is undergoing a critical re- 
examination - due in large part to its confronta- 
tion with measurements of late-type galaxy mass 
distributions indicating that dark matter is less 
concentrated than expected. In this work we 
use a recent Chandra X-ray Observatory obser- 
vation of NGC 4636 to constrain its dark matter 
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distribution, with the goal of testing whether this 
"concentration crisis" applies to this (early-type) 
elliptical galaxy. 

1.1. Chandra and Dark Matter Estimation 

Hot, extended distributions of X-ray emitting 
hot gas in hydrostatic equilibrium are crucial 
tracers of the gravitational potential in individ- 
ual elliptical galaxies gravitational lensing is 
useful only in a statistical sense for a large sam- 
ple of galaxies. The physical properties of the 
hot gas can be accurately measured out to large 
radii where dark matter dominates and optical 
techniques are infeasible. The effective resolution 
of the mass estimate essentially corresponds to 
that of the gas temperature profile (see below). 
The unprecedented Chandra angular resolution 
enables one to do X-ray imaging spectroscopy of 
the hot ISM on scales comparable to optical stud- 
ies of the stars for the first time, thus yielding the 
relative distributions of luminous and dark mat- 
ter. 

1.2. Modeling Approach 

Given a mass distribution, and the hot gas den- 
sity distribution obtained from the measured X- 
ray surface brightness distribution, we solve the 
equation of hydrostatic equilibrium to derive the 
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corresponding gas temperature profile. The to- 
tal mass model consists of three components - 
an 8 X 10^ Mq central supermassive black hole 

stars distributed as measured from Hubble 
Space Telescope and ground-based photometry, 
and a dark matter density distribution parame- 
terized by an asymptotic slope at the origin (C), 
a scale radius characterizing the transition to an 

decline (as determined by the results of nu- 
merical simulations) at large radius, and a nor- 
malization. For a given (, the dark matter scale 
and normalization, and the stellar mass-to-light 
ratio, are varied to obtain a match with the com- 
posite temperature profile measured by Chandra 
and the XMM- Newton Observatory (Figure 1). 
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Figure 1. Chandra (filled circles) and XMM- 
Newton (filled squares) temperatures, and pro- 
files for overall best-fit model (which has an r~^ 
dark matter density cusp; solid curve), and best- 
fit model with dark matter density core (broken 
curve). 

2. Results 

Models with constant mass-to-light ratios and 
those where the dark matter density inner slope, 
C, is as steep as the singular isothermal sphere 
value of —2 are clearly ruled out. We derive 



accurate constraints on the total mass distribu- 
tion from 0.7-35 kpc that are robust to the as- 
sumed value of C- The total mass increases as r^'^ 
(more quickly than the stars) to a good approxi- 
mation over this range in radii, attaining a total of 
1.5 X 10^^ Mq (corresponding to Mtot/Lv = 40 
in solar units) at the outermost point we consider 
(Figure 2). 
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Figure 2. Mass-to-light ratio, in solar units, for 
the models corresponding to the temperature pro- 
files of Figure 1 . 

There are degeneracies with respect to the rel- 
ative distributions of dark and luminous mat- 
ter. We find acceptable models with dark mat- 
ter cusps (e.g., C = 1; 1-5) and with cores {C, = 0). 
The former imply lower stellar mass to light ratios 
(the allowed range: Mstars/Ly — 2.4 — 6.6, consis- 
tent with studies of the stellar population) and, in 
some cases, the predominance of dark matter all 
the way into the nucleus of the galaxy. The dark 
matter fraction within the half-light radius (~ 8 
kpc) can be as high as 0.8 and "central" (actually, 
the average over the inner 700 pc) dark matter 
density as high as 4.0 M© pc~^; the cored mod- 
els define lower limits to these quantities of 0.5 
and 0.16 M© pc~^, respectively. Even the lowest 
allowed central density is an order of magnitude 
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greater than the previously estimated values for 
other types of galaxies (Figure 3) that precipi- 
tated the consideration of alternatives to CDM. 
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Figure 3. Average enclosed dark matter den- 
sity distribution corresponding to the tempera- 
ture profiles of Figure 1. The typical value esti- 
mated for spiral (and other) galaxies is indicated. 



3. Implications for the Nature of Dark 
Matter 

Comparison of non-parametric measures of 
the dark matter concentration with expectations 
based on standard CDM numerical simulations 
indicate that dark matter in NGC 4636 is more 
compact than average for its mass scale, but 
within the expected scatter. This is consistent 
with the findings of statistical studies of the grav- 
itational lensing properties of elliptical galaxies, 
though it appears that NGC 4636 is unusually 
dark matter dominated. 

Many of the proposed scenarios for reducing 
the dark matter concentration in late-type galax- 
ies are at odds with our results. It is instruc- 
tive to consider the two classes of models illus- 
trated in the above Figures from this perspective 
- those with flat dark matter cores, and those 



with cuspy dark matter cores (that provide bet- 
ter fits to the X-ray data) interpreted as contract- 
ing from initially flat cores due to baryonic in- 
fall. For the flat-core models, the high central 
dark matter mass density and large cores are con- 
trary to the expected scaling relations for self- 
interacting dark matter (SIDM) and other mod- 
els where dark matter structure is driven by dark 
matter particle interaction. For the cuspy-core 
models the central phase-space density, conserved 
during adiabatic contraction, is too high. Alter- 
natively, the velocity-dependent cross section can 
be fine-tuned so that interactions are ineffective 
at the higher mass (and velocity) scale of ellipti- 
cals; although here too, one may run into phase 
space difficulties if ellipticals form from mergers 
of lower mass systems with interaction-induced 
dark matter cores. 

Scenarios where the mechanism for reducing 
the dark matter concentration becomes less ef- 
fective with increasing mass scale are not ruled 
out, since ellipticals such as NGC 4636 represent 
the most massive galaxies with the deepest po- 
tential wells. These include models that invoke 
warm dark matter, or the expansion of dark mat- 
ter due to coupling with a powerful protogalactic 
starburst-driven baryonic outflow. If this transi- 
tional mass is indeed on the giant elliptical galaxy 
scale, cuspy dark matter distributions in galaxy 
clusters are implied. 

Our result that the central dark matter den- 
sity in NGC 4636 is > 8 - 200 x the typical value 
for less massive galaxies contradicts the naive ex- 
pectations of bottom-up hierarchical clustering 
where the most massive systems form latest and 
reflect the relatively low average density in the 
universe at that epoch. A cold dark matter ini- 
tial perturbation spectrum that is tilted or oth- 
erwise lacking in small-scale power may help ex- 
plain Ly-a forest data Q, but likely exacerbates 
this contradiction. 

Recent work on the relation between the merg- 
ing histories of galaxies and their morphology |g] 
may be the key to resolving the dark matter con- 
centration problem. A significant dispersion ~ in 
dark matter concentration is predicted by numer- 
ical simulations - at any given mass range, but 
particularly on galaxy scales. This results from 
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the stochasticity of the merging process, and in- 
troduces the following bias. Galaxies of low cen- 
tral dark matter density represent relatively re- 
cently formed systems, and/or particularly fragile 
galaxies that experienced relatively tranquil as- 
sembly histories these may correspond to disk 
galaxies where low dark matter central densities 
are indeed measured. Conversely, giant elliptical 
galaxies such as, or perhaps particularly, NGC 
4636 may form at relatively high redshifts and un- 
dergo exceptionally prominent merger histories. 

4. Concluding Remarks 

Possible observational inconsistencies of the 
(non-interacting) CDM paradigm resulted in a re- 
newed scrutiny and a proliferation of alternative 
models. Careful analysis of the observational sit- 
uation in light of more detailed, and more deeply 
understood, theoretical models, as well as the 
shortcomings of the proposed alternatives seem 
to indicate that CDM is withstanding these re- 
cent challenges to remain the most viable model 
of structure in the universe (see J. Primack's con- 
tribution to this volume) . Our results on the high 
concentration of dark matter in NGC 4636 evi- 
dently strengthens the case for CDM, while in- 
troducing other puzzles - such as why NGC 4636 
seems to be exceptionally dark matter dominated. 
Although one must exercise caution in general- 
izing from our results, the most stringent con- 
straints on alternative models for dark halo struc- 
ture, that presumably universally apply, emerge 
from studying such a, possibly extreme, system. 
Results, in progress, of similar investigations for 
other galaxies (e.g., NGC 1399, NGC 4472) are 
anticipated with great interest. 
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